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摘 要 : 利用 某 一 风电 场 区 域内 33 人 台风 机 观测 资料 以 及 欧洲 中 期 天 气 预报 中 心 (ECMWF) 的 100 m 
高 度 风 速 预报 产品 ,引进 最 优 训 练 期 方案 ,通过 一 元 线性 回归 方法 先 对 模式 的 100m 高 度 风速 预报 
进行 订正 ,再 分 别 采 用 完全 预报 方法 (PP 法 ) 和 模式 输出 统计 预报 方法 (MOS 法 ) 进 行 风 电功率 预报 


的 对 比 研 究 。 结 


表明 :ECMWF 的 100m 高 度 风速 预报 值 与 风机 轮 载 高 度 处 风速 观测 值 之 间 的 偏 


差 并 不 大 ,订正 后 的 预报 误差 进一步 减 小 。 将 订正 后 的 风速 预报 代入 PP 法 和 MOS 法 建立 的 风电 
功率 预报 方程 ,可 以 显著 减 小 PP 法 的 预报 误差 ,但 是 并 没有 改进 MOS 法 的 预报 结果 。 与 利用 风速 
预报 订正 产品 进行 风电 功率 预报 的 PP 法 相 比 ,MOS 法 可 以 省 去 风速 预报 的 订正 环节 而 略 胜 一 筹 ， 
简化 了 业务 流程 。 这 2 种 方法 计算 的 33 合 风机 站 3~72h 预 报 的 均 方 根 误差 和 平均 绝对 误差 分 别 


在 17%~25% 和 和 11%~18% 之 间 , 具 有 业务 应 用 价值 。 
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在 全 球 气候 变 暧 和 能 源 日 趋 紧张 的 背景 下 , 风 
力 发 电 已 成 为 世界 上 最 引 人 注 目的 新 型 能 源 ” , 然 
而 风能 资源 存在 较 大 的 间 欣 性 和 波动 性 ,大 规模 的 
风电 接 入 会 影响 电力 系统 的 安全 和 稳定 运行 A 
此 ,准确 的 风速 和 风电 功率 预测 是 保障 电能 质量 和 
风电 合理 调配 的 迫切 需求 ” 。 

风电 功率 与 风速 的 三 次 方 成 正比 “ ,对 不 同 区 
域 风速 的 变化 规律 进行 分 析 ”” ,对 做 好 风能 资源 评 
估 具 有 重要 意义 ,而 风电 功率 预测 系统 的 核心 部 分 
就 是 对 风 的 预报 。 数 值 天 气 预报 是 风电 场 风 速 预 
报 的 主要 文 撑 手 段 ,美国 大 气 研 究 中 心 和 深 州 大 学 
人 研发 的 MM5(Mesoscale model) 丹麦 的 KAMM(Karl- 
sruhe atmospheric mesoscale model) , 美国 国家 海洋 
和 大 气管 理 局 的 RAMS(Regional atmospheric model- 
ling system) 和 美国 国家 大 气 人 研究 中 心 的 WRF (The 
weather research and forecasting mode]) 等 中 尺度 数值 
模式 在 风电 场 风速 预报 中 发 挥 了 重要 作用 ,但 
是 很 多 省 ,市 级 气象 部 门 不 具备 直接 获取 中 太 度 数 
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值 模式 产品 的 能 力 ,给 电力 气象 服务 工作 带 来 了 较 
大 困难 。 

ECMWF 是 气象 部 门 常 用 的 具有 较 高 预报 水 平 
的 全 球 数值 模式 "", 可 以 提供 100 m 高 度 (比较 接近 
风机 轮 载 高 度 ) 的 风速 预报 ,然而 目前 利用 该 模式 
产品 进行 风电 场 风 速 和 功率 预测 的 研究 还 比较 少 
见 。 在 数值 天 气 预 报 中 , 近 地 面 风 的 预报 难度 远 远 
高 于 具有 高 连续 性 变化 的 温度 的 预报 ,这 主要 因为 
近 地 面 风 不 满足 地 转 平衡 , 除 受 大 尺度 的 系统 性 天 
气 影响 外 ,还 受 地 形 、 下 垫 面 热 动 力作 用 以 及 边界 
层 中 夹 卷 过 程 的 影响 ,具有 较 强 的 局 地 性 、 瞬 变性 
和 间 向 性 ,因此 ,无 论 何 等 分 辩 率 的 数值 模式 产 
品 都 很 难 精确 到 具体 描述 风机 风 场 的 变化 。 数 值 
模式 预报 结合 统计 订正 是 提高 风速 预报 准确 率 的 
常用 方法 ,如 人 工 神经 网 络 支持 向 量 机 "7 、 最 
小 二 乘法 卡尔 曼 滤波 ”以 及 相似 误差 订正 
等 。 一 元 线性 回归 方法 作为 最 常用 的 统计 方法 之 
一 ,近年 来 被 气象 部 门 广泛 应 用 于 气温 预报 业务 ””， 
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但 是 在 风速 预报 业务 中 的 应 用 并 不 多 见 。 利 用 统 
计 方 法 建立 预报 方程 时 ,可 以 采用 最 优 训练 期 来 提 
高 预报 质量 , 当 训 练 期 短 于 最 优 训 练 期 长 度 时 , 预 
报 误差 并 未 达到 最 小 ,而 超过 此 训练 期 长 度 后 , 预 
报 误差 并 没有 很 大 的 降低 ,反而 由 于 训练 期 的 增长 
增加 了 计算 量 和 计算 时 间 '”。 吴 启 树 等 ”利用 最 
优 训 练 期 方案 ,显著 提高 了 一 元 线性 回归 方法 对 模 
式 气温 预报 的 订正 能 力 。 但 是 一 元 线性 回归 方法 
也 有 局 限 性 ,一 方面 ,如 果 预 报 对 象 和 预报 因子 之 
间 的 相关 性 不 显著 , 则 统计 预报 结果 不 会 太 理 想 ， 
另 一 方面 ,统计 预报 结果 的 精细 化 水 平 决 定 于 预报 
因子 的 数据 分 辨 率 。 

完全 预报 方法 (PP 法 ) 和 模式 输出 统计 预报 方 
法 (MOS 法 ) 是 2 种 比较 常用 的 数值 预报 产品 使 用 方 
法 。PP 法 利用 历史 资料 中 预报 因子 的 观测 资料 与 
同时 间 的 预报 对 象 建立 统计 关系 ,而 MOS 法 利用 历 
史 资 料 中 预报 因子 的 数值 预报 产品 ,与 预报 时 效 对 
应 时 刻 的 预报 对 象 建立 统计 关系 ,在 预报 时 ,这 两 
种 方法 都 将 预报 因子 的 数值 预报 结果 代入 所 建立 
的 预报 方程 做 出 预报 。PP 法 的 预报 结果 除了 包含 
建立 统计 方程 时 的 拟 合 误差 外 ,还 包含 数值 预报 模 
式 的 误差 ,而 MOS 法 将 数值 预报 模式 的 误差 及 局 地 
气候 影响 等 因素 直接 统一 到 建立 预报 方程 的 拟 合 
误差 中 ,从 而 减少 了 误差 源 下 。 如 果 先 对 模式 的 风 
速 预报 进行 订正 ,再 代入 PP 法 和 MOS 法 建立 的 风 
电功率 预测 模型 ,无 疑 可 以 减 小 PP 法 的 风电 功率 预 
报 误差 ,但 是 否 一 定 会 减 小 MOS 法 的 预报 误差 呢 ? 
另外 ,利用 MOS 法 建立 的 风电 功率 预报 方程 是 否 一 
定 优 于 PP 法 呢 ? 

本 文 利 用 ECMWF 的 100 m 高 度 风速 预报 ,插值 
得 到 风机 位 置 的 风 场 预报 ,然后 采用 基于 最 优 训练 
期 的 一 元 线性 回归 方法 订正 风速 预报 场 , 最 后 通过 
PP 法 和 MOS 法 进行 风电 功率 预报 的 对 比 研 究 。 
ECMWF 的 100 m 高 度 风 速 预报 是 逐 3h 或 者 6h 间 
隔 的 ,与 电力 生产 调度 上 逐 15 min 间隔 的 预报 要 求 
有 较 大 差距 ,但 是 可 以 用 来 预测 风电 场 的 风速 趋 
势 ,并 且 本 研究 旨 在 提供 一 种 将 模式 的 风速 预报 合 
理应 用 到 风能 预报 业务 的 有 效 方法 ,如果 以 更 高 分 
辩 率 的 模式 产品 作为 研究 资料 , 则 可 以 得 到 更 高 分 
辩 率 的 风速 和 风电 功率 预测 产品 。 


1 资料 与 处 理 


观测 资料 来 自 2016 年 1 月 1 日 一 2018 年 12 月 


31 日 陕西 省 某 一 风电 场 区 域内 1 台 测 风 塔 的 80 m 
风向 观测 和 33 台 风机 的 风速 有功 功率 观测 ,风机 
轮 坑 高 度 为 80 m, 额 定 功率 为 1500 kW ,额定 风速 为 
10 ms ,切入 风速 为 3 ms , 切 出 风速 为 20 m-s, BK 
据 采 集 频 率 为 5 min。 为 了 保证 研究 结果 的 可 靠 性 ， 
剔除 风机 观测 资料 中 的 明显 不 合理 值 '” ,主要 包括 
剔除 风速 大 于 5 ms 而 有 功 功 率 小 于 50 kW 的 时 
刻 的 观测 值 ; 剔 除 有 功 功率 大 于 额定 功率 1500 kW 
时 刻 的 观测 值 。 

模式 资料 采用 2016 年 1 月 1 日 一 2018 年 12 月 
31 H ECMWF 的 100 m 高 度 上 东西 风 (V) 和 南北 风 
(四 的 预报 。 由 于 ECMWF 数 据 传 输 至 业务 系统 的 时 
间 比 起 报时 间 滞 后 6h 以 上 ,因此 ,采用 比 起 报时 间 
提前 12h 的 预报 , 即 采 用 前 一 日 20:00( 当 日 08:00) 
起 报 的 15~240 bh 的 风速 ,制作 当日 08:00(20:00) 起 
报 的 3~228ph 的 风机 轮 载 高 度 处 的 风速 和 风电 功 
率 ,预报 时效 为 3~60h 逐 3h 间 隔 、66~228h 逐 6h 
间隔 ,水 平分 辨 率 为 0.125°x0.125°*。 由 于 分 辨 率 不 
够 精细 ,ECMWF 模式 不 适用 于 复杂 地 形 ,而 该 风电 
场 的 场 址 区 海拔 高 度 约 1400 m, Jy E EIE A E Be e 
Ih ,场地 较为 开阔 ,地 势 起 伏 不 大 (图 1) ,为 本 文采 
H ECMWF 的 100 m 高 度 风速 预报 作为 研究 资料 提 
供 了 依据 。 


N 
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图 例 
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e 风机 站 位 置 
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高 : 1591 
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图 1 风电 场 及 其 周边 地 区 的 海拔 高 度 .ECMWF 的 
格 点 和 风机 站 的 位 置 分 布 示 意图 
Fig. 1 Altitude of the wind farm and its surround area, 
and location distributions of wind turbines and lattice 
points of the ECMWF mode 
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为 了 使 观测 资料 与 模式 资料 相 匹 配 ,对 研究 资 
料 进行 如 下 处 理 :(1) 利用 模式 的 100m 高 度 上 东西 


UCU) AUPE UV) AIMAR HI AK WS = NU’ +V, 


计算 得 到 100 m 高 度 全 风速 (WS, 单 位 ms ') 预 报 ; 
(2) 采用 双 线 性 插值 方法 “将 模式 的 格 点 值 插值 到 
站 点 上 , 即 , 每 个 风机 站 的 值 都 通过 周围 4 个 格 点 值 
在 纬 问 和 经 和 癌 方 向 分 别 进行 一 元 一 次 线性 插值 后 
得 到 ,模式 的 格 点 和 33 台 风机 的 位 置 分 布 如 图 1 所 
示 ;(3) 将 观测 资料 处 理 为 与 模式 的 预报 时 效 相 匹 
配 的 定时 资料 ,例如 ,2017 年 1 月 1 日 08:00 起 报 的 
3h、6h、12h 等 预报 时 效 的 风速 预报 ,与 2017 年 1 月 
1 日 11:00、14:00 、17:00 等 整 点 时 刻 的 风速 观测 值 相 
对 应 ,其 它 预报 时 效 依 此 类 推 。 


2 研究 方法 


2.1 一 元 线性 回归 方法 

利用 模式 预报 和 风电 场 观 测 资料 ,计算 得 到 模 
式 的 风速 预报 值 与 风机 的 风速 和 输出 功率 观测 值 
的 相关 系数 分 别 为 0.47 和 0.48 ,通过 了 95% 的 显著 
性 检验 ,为 利用 一 元 线性 回归 方法 进行 风电 场 的 风 
速 和 输出 功率 预测 提供 了 理论 依据 。 由 于 预报 服 
务 人 员 开 始 制作 预报 时 ,预报 日 前 1 天 的 实况 资料 
不 全 ,因此 ,训练 样本 采用 预报 日 前 1 天 之 前 的 Nd 
(训练 期 长 度 ) 资 料 , 随 预报 日 滑动 取样 。 在 训练 
期 ,分 别针 对 每 台风 机 站 建立 各 预报 时 次 的 风机 轮 
载 高 度 处 的 风速 或 者 风电 功率 预报 的 一 元 线性 回 
归 方 程 : 


y=axt+b (1) 
式 中 :7 为 预报 对 象 ;x 为 预报 因子 ; 当 训 练 样本 数 >4 
时 ,通过 最 小 二 乘法 计算 斜率 (a) 和 截 距 (2)。 
2.2 最 优 训 练 期 及 其 评估 周期 
以 10~100 d 逐 日 作为 训练 期 日 数 (N) ,以 一 段 
时 期 作为 最 优 训 练 期 的 评估 周期 ,根据 该 评估 周期 
内 33 台风 机 的 3~228h 风 速 (或 者 风电 功率 ) 预 报 的 
均 方 根 误差 随 训 练 期 日 数 的 变化 ,确定 最 优 训练 期 
日 数 (N')。 由 于 需要 100 d( 最 大 的 训练 期 日 数 ) 的 前 
置 历史 资料 作为 训练 期 ,因此 ,对 于 风速 预报 ,2016 
年 1 月 1 日 一 4 月 10 日 的 数据 只 能 作为 训练 期 数据 ， 
实际 的 风速 预报 从 2016 年 4 月 11 日 开始 ,风电 功率 
预报 是 在 风速 预报 的 基础 上 开展 的 ,对 于 风电 功率 
预报 ,2016 年 4 月 11 日 一 7 月 19 日 的 数据 只 能 作为 训 
练 期 数据 ,实际 的 功率 预报 从 2016 年 7 月 20 日 开始 。 


最 优 训练 期 日 数 N’' 应 该 是 预报 日 前 给 定 的 。 
本 文 设 计 了 5 种 在 预报 日 前 利用 不 同 评估 周期 滑动 
确定 最 优 训练 期 日 数 的 方案 ( 表 1) ,由 于 需要 365 d 
(最 大 的 评估 周期 ) 的 前 置 历史 资料 作为 最 优 训练 
期 日 数 的 评估 周期 ,风速 预报 和 风电 功率 预报 的 检 
验 分 别 从 2017 年 4 月 11 日 和 7 月 20 日 开始 。 将 2017 
年 7 月 20 日 之 后 的 数据 分 为 2017 年 7 月 20 日 一 
2018 年 2 月 28 日 和 2018 年 3 月 1 日 一 12 月 31 日 两 
部 分 ,前 者 用 来 对 方案 1~5 的 风速 (风电 功率 ) 预 报 
效果 进行 评估 ,从 而 选 出 最 佳 的 评估 周期 方案 ,后 
者 用 来 对 最 佳 评估 周期 方案 下 的 风速 (风电 功率 ) 
预报 效果 进行 检验 。 


R1 最 优 训练 期 的 评估 周期 方案 
Tab.1 Evaluation cycle schemes of the 
optimal training period 


方案 评估 周期 

预报 日 前 1 天 之 前 的 365 d(1 年 ) 
预报 日 前 1 天 之 前 的 182 d(6 个 月 ) 
预报 日 前 1 天 之 前 的 91 d(3 个 月 ) 
预报 日 前 1 天 之 前 的 30 d(1 个 月 ) 
预报 日 前 1 天 之 前 的 10d 


nan Ae ù N =e 


2.3 风速 预报 订正 方法 和 检验 指标 

逐 预 报时 次 滑动 取 预 报 日 前 1 天 之 前 最 优 训练 
期 ( W' ) 的 风速 预报 和 观测 资料 作为 训练 样本 ,利用 
一 元 线性 回归 方法 ,分 别针 对 每 台风 机 站 建立 风速 
预报 订正 方程 ,将 模式 的 100 m 高 度 风 速 预报 订正 
到 风机 轮 载 高 度 上 。 当 风速 预报 的 订正 值 <0 时 , 令 订 
正 值 为 0。 为 了 客观 评估 本 文 对 模式 100 m 高 度 风 速 
预报 的 订正 效果 ,风速 预报 的 检验 指标 不 仅 包括 平 
均 绝 对 误差 (War) 均 方 根 误差 (fwss) 和 订正 结果 相 
对 模式 改善 百分比 (PM, 以 下 简称 风速 改善 百分比 )， 
还 包括 系统 误差 (SYS_E) 和 随机 误差 (RAN_E)'”,， 
计算 公式 如 下 : 


SIRA (2) 
i=l 
1 f 2 
Vane = D -y,) (3) 
i=l 
Vive =V 
PM = 7 x 100% (4) 


MAE, 


SYS_E= | Dg 5 ya)| (5) 


i= 


RAN_E= ee E y,) ~ D ~ ya) (6) 
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GLP sy, y, 分 别 指 第 ;次 预报 的 预报 值 和 观测 值 ， 
预报 对 象 为 风速 ;n 为 参与 检验 的 总 预报 次 数 ; 
os > Vue, 分 别 指 风速 预报 订正 前 、 后 的 平均 绝对 
误差 ;PM>0(PM<0) 表 明 订 正 效果 为 正 ( 负 )。 
2.4 风电 功率 预报 方法 和 检验 指标 

逐 预 报时 次 滑动 取 预 报 日 前 1 天 之 前 最 优 训练 
期 ( NW' ) 的 风机 功率 观测 和 风速 观测 (或 者 风速 预 
报 ) 资 料 作为 训练 样本 ,基于 一 元 线性 回归 方法 , 利 
用 PP 法 (或 者 MOS 法 ) 建 立 每 一 台风 机 站 的 功率 预 
测 方程 ,将 订正 前 、 后 的 风速 预报 分 别 代 和 人 预报 方 
程 制作 风电 功率 预报 ,预报 产品 分 别 记 作 PP1、PP2 
(ak MOSI.MOS2), B4«Q/GDW10588—2015 风 
电功率 预测 功能 规范 》, 风 电功率 预报 的 检验 指标 
包括 平均 绝对 误差 (Puws) 和 均 方 根 误差 (Pawss) , 计 
算 公 式 如 下 : 


= 
N 
© 


(a) 方案 1 
V 
R 90 
K 
5 60 
Kay 
I 30 
Y GF QA SS VY yyy Q 
Tà a a Vv NO N a a 
YY YY YY YY LY LY D 
FP TT HP HP TT 7 5 
日 期 (年 -月 -日 ) 
120 
(c) 方案 3 
7 
R 
KN 
= 60 f 
as 
R 30 
Y 8 Q SS VY ££ YY Lf 
~ a a Vv NYO N m wa 
Yr TF TY YT LY LY © 
小 路 小 小 小 小 中 F 
日 期 (年 -月 -日 ) 
120 
mj 
= 90 
$ 
= 60 
Ka 
i 30 
XY CS SF GT NX KY LY 
a > “> Vv NO N wa a 
YY Sy YY TY LY LY 
rr PK FFY 
日 期 (年 -月 -日 ) 


te Vo, 
Py = | (7) 


(8) 


AP: y,、y, 分 别 指 第 i 次 预报 的 预报 值 和 观测 值 ， 
为 参与 检验 的 总 预报 次 数 。 


3 风速 预报 订正 结果 


3.1 最 优 训 练 期 的 确定 

利用 最 优 训练 期 的 评估 周期 方案 1~5 ,分别 将 
2017 年 7 月 20 日 一 2018 年 2 月 28 日 模式 的 100m 高 
度 风 速 预 报 订 正 到 风机 轮 载 高 度 上 。 从 最 优 训练 
期 日 数 来 看 (图 2) ,方案 5 的 稳定 性 最 差 ,这 与 较 短 


me 
= 90 
€ 
5 60 
aS 
IR 30 
© X er DYP 
N S Q Y SF YJ y 
Ma a ae Vv YO N a af 
yy Fy TY LY LY TY LY eD 
PF KP HP HP HP KF F 
日 期 (年 -月 -日 ) 
120| (a) WH 
= 90 
= 
S 
= 60 
= 
IE 30 
~ E S SS LS y 
“> a ae NYNO N N wa wa 
“S Fy TY TY LY YY &\ Sy 
PP TT HP TT TT F 个 
日 期 (年 -月 -日 ) 


图 2 2017 年 7 月 20 日 一 2018 年 2 月 28 日 风速 预报 订正 的 逐日 最 优 训 练 期 日 数 
Fig. 2 Daily optimal training days for correcting wind speed forecast from 20" July 2017 to 28" February 2018 
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周期 确定 的 评估 周期 容易 受 转折 性 天 气 过 程 的 影 
响 有 关 汪 。 但 是 从 均 方 根 误 差 来 看 (图 3) ,方案 5 订正 
后 的 风速 预报 的 均 方 根 误差 最 小 。 另 外 ,本 文 还 参考 
文献 [23] 中 的 最 优 训练 期 方案 ,分 别 以 1 年 .6 个 月 、 
3 个 月 和 1 个 月 作为 评估 周期 (M) ,应 用 上 1 年 逐 朵 
的 预报 评估 确定 该 年 逐 导 的 最 优 训 练 期 日 数 ,对 模 
式 的 风速 预报 进行 订正 ,并 与 方案 5 的 风速 预报 作 
比较 ,研究 结果 再 次 表明 方案 5 的 预报 误差 最 小 ,与 


4.5 


(a) 订正 前 
8 4.0 
a 
mK 3.5 
T 
p 3.0 
2 


15 30 45 60 90 120 150 180 210 
预报 时 效 /h 


文献 [23] 研 究 认 为 以 较 短 周期 为 评估 周期 滑动 确定 
的 最 优 训练 期 不 能 进一步 提高 气温 预报 效果 的 结 
论 相 反 , 这 可 能 因为 本 文 的 订正 对 象 为 风速 预报 ， 
而 气温 和 风速 是 2 种 变化 特征 截然 不 同 的 气象 要 素 。 
3.2 预报 误差 分 析 

利用 方案 5 逐日 滑动 确定 最 优 训 练 期 日 数 ,对 
2018 年 3 月 1 日 一 12 月 31 日 风电 场 33 台 风机 站 的 
风速 进行 预报 并 检验 ,图 4 是 08:00 起 报 的 33 台风 


(b) 订正 后 一 方案 1 — 方案 4 
-一 方案 2 — 方案 5 
-一 方案 3 


均 方 根 误差 /ms 


15 30 45 60 90 120 150 180 210 
预报 时 效 /h 


图 3 2017 年 7 月 20 H—20184F 2A 28 日 风速 预报 订正 前 、 后 的 均 方 根 误差 


Fig.3 Root mean square error of wind speed forecast before and after correction from 20" July 2017 to 28" February 2018 
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图 4 2018 年 3 月 1 日 一 12 月 31 日 08:00 起 报 的 33 台风 机 站 平均 的 风速 预报 值 与 观测 值 的 对 比 及 预报 误差 分 析 


Fig.4 Comparison between observation and prediction of wind speed from 08:00 BST and the 


analysis of its forecast error at 33 wind turbines from 1“ March to 31" December 2018 
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机 站 平均 风速 预报 的 检验 结果 (20:00 起 报 的 图 
略 )。 从 图 4a~b 来 看 ,模式 的 100 m 高 度 风 速 预报 
值 虽 然 比 风机 轮 载 高 度 处 的 风速 观测 值 偏 大 一 些 ， 
但 是 却 可 以 较 好 地 预报 出 风速 的 日 变化 和 逐 月 变 
化 趋势 ,订正 后 的 风速 预报 误差 显著 减 小 ,预报 值 
和 观测 值 的 变化 曲线 几乎 重合 。 从 图 4c~f 来 看 ,3~ 
228h 预 报时 效 ,模式 的 风速 预报 的 系统 误差 和 随机 
误差 分 别 从 订正 前 的 0.9~1.6 ms* 和 1.9~3.7 mes BE 
低 为 订正 后 的 0.0~0.1 ms #0 1.6~2.6 m-s”, IIJ AR 
误差 和 平均 绝对 误差 分 别 从 订正 前 的 2.2~3.9 mes! 
和 1.7~3.2 ms 降低 为 订正 后 的 1.6~2.6 ms il 
1.2~2.1 ms 。 从 图 4g 来 看 ,订正 前 ,风速 预报 的 误 
差 主 要 出 现在 -2.0~0.0 m-s',0.0~2.0 ms 和 2.0~ 
4.0 ms ,出 现 概率 分 别 为 21%、 31% 和 24% , 正 误差 
(预报 值 大 于 观测 值 ) 比 例 为 69% ;订正 后 ,风速 预 
报 误 差 主 要 出 现在 -2.0~0.0 ms"' 和 0.0~2.0 mss", 
出 现 概率 分 别 为 32% 和 37%, 正 误差 比例 为 53%, 订 
正 后 的 预报 误差 和 正 、 负 误差 比例 之 间 的 差距 都 小 
于 订正 前 。 从 图 4h 来 看 ,订正 后 的 风速 改善 百分比 
在 24%~41% 之 间 ,1 日 中 ,风速 改善 百分比 的 极 大 
值 和 极 小 值 分 别 出 现在 中 午 和 凌晨 ,与 风速 预报 值 
与 观测 值 的 相关 系数 的 极 大 值 和 极 小 值 的 出 现时 
段 一 致 。 网 $ 进 一 步 检 验 了 逐 人 台风 机 站 的 风速 预报 
订正 效果 ,结果 表明 ,33 台风 机 站 的 风速 预报 误差 
均 得 到 了 明显 改善 ,3~72h 预 报 的 平均 绝对 误差 从 
订正 前 的 1.8~2.4 ms 降低 为 订正 后 的 1.2~1.6 m+s”, 
风速 改善 百分比 在 12%~50% 之 间 , 在 风速 预报 误 
差 越 大 的 风机 站 改善 效果 越 好 。 下 面 将 采用 方案 5 
逐日 滑动 确定 最 优 训练 期 日 数 ,对 2016 年 4 月 11 
日 一 2018 年 12 月 31 日 风电 场 33 台 风机 站 上 的 模式 
风速 预报 进行 订正 ,然后 将 订正 结果 作为 风电 功率 
预报 的 研究 资料 。 


4 风电 功率 预报 结果 分 析 


4.1 最 优 训 练 期 的 确定 

分 别 利 用 最 优 训 练 期 的 评估 周期 方案 1~5, 采 
用 PP 法 和 MOS 法 制作 2017 年 7 月 20 日 一 2018 年 2 
月 28 日 的 风电 功率 预报 产品 PP1、PP2 和 MOS1、 
MOS2, 并 检验 最 优 训练 期 日 数 和 均 方 根 误差 。 以 
MOS1 为 例 ( 图 6) ,方案 1~5 逐 日 滑动 确定 的 最 优 训 
练 期 日 数 的 不 稳定 性 随 着 评估 周期 的 变 短 而 增强 ， 
但 是 方案 1~5 的 均 方 根 误差 几乎 重合 。 可 见 ,以 较 


£27) @ 订正 前 的 风速 预报 


3 6 9 12 15 18 21 24 30 33 
站 号 序列 


平均 绝对 误差 


3 6 9 12 15 18 21 24 30 33 
站 号 序列 


(c) 风速 改善 百分比 


风速 改善 百分比 /% 


3 6 9 12 15 18 21 24 30 33 
站 号 序列 


图 5 2018 年 3 月 1 日 一 12 月 31 日 08:00 起 报 的 
各 风机 站 3~72h 风 速 预报 订正 前 、 后 的 
平均 绝对 误差 和 风速 改善 百分比 
Fig. 5 Mean absolute error of wind speed prediction before 
and after correction and the percentage of improvement in root 
mean square error of wind speed forecast after correction 
during 3 to 72 forecast hours form 08:00 BST in different 
wind turbines from 1“ March to 31" December 2018 


长 周期 为 评估 周期 滑动 确定 的 最 优 训练 期 的 稳定 
性 虽然 优 于 较 短 周期 ,但 是 并 不 能 进一步 提高 预报 
效果 。 针 对 PP1、PP2 和 MOS2 做 如 上 分 析 , 所 得 结 
论 与 MOS1 一 致 。 由 于 方案 5 需要 的 训练 样本 最 少 ， 
有 利于 业务 化 应 用 ,所 以 将 方案 5 作为 最 优 训 练 期 
的 最 佳 评 估 周 期 方案 ,制作 2018 年 3 月 1 日 一 12 月 
31 日 各 风机 站 的 风电 功率 预报 产品 PP1、PP2 和 
MOS1 MOS2 。 
4.2 预报 误差 分 析 

图 7 对 比 检验 了 2018 年 3 月 1 日 一 12 月 31 日 
08:00 起 报 的 PP1、PP2、MOS1 和 MOS2 的 平均 功率 
预报 (20:00 起 报 的 图 略 ) ,该 平均 功率 预报 值 通过 将 
33 台风 机 在 同一 预报 时 效 的 功率 相 加 后 除 以 风机 
台数 得 到 。 从 图 7a~b 来 看 ,PP1、PP2、MOS1 和 
MOS2 都 可 以 较 好 地 预报 出 风电 场 功 率 输出 的 日 变 
化 和 逐 月 变化 趋势 ,但 是 PP1 的 预报 值 明 显 大 于 观 
测 值 ,这 与 订正 前 的 模式 100 m 高 度 风 速 预 报 值 大 
于 风机 轮 载 高 度 处 的 风速 观测 值 有 关 ,而 PP2、 
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图 6 2017 年 7 月 20 日 一 2018 年 2 月 28 日 方案 1~5 的 风电 功率 预报 的 逐日 最 优 训练 期 日 数 和 均 方 根 误差 
Fig.6 Daily optimal training days and the root mean square error of wind power prediction 
of schemes 1-5 from 20" July 2017 to 28" February 2018 


MOS1 和 MOS2 的 预报 曲线 与 观测 曲线 几乎 重合 。 
从 图 7c~e 来 看 ,3~228 h 预 报时 效 ,PP2、MOS1 和 
MOS2 的 均 方 根 误差 均 为 16%~31% ,平均 绝对 误差 
依次 为 10%~25% 、10%~25% 和 10%~26% ,绝对 误差 
三 20% 的 概率 依次 为 68% 、68% 和 67%。 男 外 ,这 4 
种 风电 功率 预报 产品 的 误差 极 大 值 ( 极 小 值 ) 都 出 
现在 上 午 ( 凌 晨 ) ,与 风速 预报 误差 的 极 大 值 ( 极 小 
值 ) 的 出 现时 段 一 致 。 从 以 上 数据 来 看 , 先 对 模式 
的 100 m 高 度 风速 预报 进行 订正 ,再 代入 PP 法 和 
MOS 法 建立 的 风电 功率 预报 方程 ,可 以 显著 减 小 PP 
法 的 预报 误差 ,但 是 并 没有 改进 MOS 法 的 预报 结 
果 。 与 利用 风速 预报 订正 产品 进行 风电 功率 预报 
的 PP 法 相 比 ,MOS 法 可 以 省 去 风速 预报 的 订正 环节 
而 略 胜 一 筹 ,简化 了 业务 流程 。 

图 8 以 08:00 起 报 的 MOSI1 为 例 ,检验 了 各 风机 


站 的 功率 预报 效果 。 从 图 8 来 看 ,MOS1 的 风电 功率 
预报 值 随 风机 站 号 的 变化 趋势 与 观测 值 一 致 ,3~72h 
预报 的 均 方 根 误差 和 平均 绝对 误差 分 别 在 17% ~ 
25% 和 11%~18% 之 间 ,预报 值 和 观测 值 人 篇 大 的 站 预 
报 误差 也 偏 大 一 些 。 对 PP2 和 MOS2 做 如 上 分 析 ， 
所 得 结论 与 MOS1 一 致 。EWind 和 Prediktor 在 美国 
同时 为 2 个 大 型 风电 场 发 布 的 风电 功率 预报 的 平均 
绝对 误差 分 别 为 14%-~22% 和 14% ~18% , Hh JI Ae 
利用 神经 网 络 方法 计算 的 整个 风电 场 72h 风 电 功 
率 预 报 的 均 方 根 误 差 在 21%~24% 之 间 ,与 以 上 研究 
相 比 ,PP2、MOS1 和 MOS2 的 预报 误差 并 不 大 ,具有 
业务 应 用 价值 。 


5 讨论 


本 文 的 风电 功率 预报 主要 考虑 了 风速 的 影响 ， 
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图 7 2018 年 3 月 1 日 一 12 月 31 日 08:00 起 报 的 33 台风 机 站 的 平均 功率 预报 值 与 观测 值 的 对 比 及 预报 误差 分 析 


Fig. 7 Comparison between observation and prediction of wind power from 08:00 BST and the analysis 


of its forecast error at 33 wind turbines from 1“ March to 31" December 2018 
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图 8 2018 年 3 月 1 日 一 12 月 31 日 08:00 起 报 的 3~72h 风 电功率 预报 的 逐 台 风机 站 检验 
Fig. 8 Quality test of wind power forecast during 3 to 72 forecast hours from 8:00 BST in 
different wind turbines from 1“ March to 31" December 2018 


而 风机 功率 输出 还 与 空气 密度 密切 相关 ,但 空气 密 
度 比 风速 对 风机 功率 输出 的 影响 小 得 多 ,并 且 该 风 
电场 82% 的 观测 时 刻 的 空气 密度 与 其 多 年 平均 值 
的 相对 偏差 小 于 5%。 因 此 本 文 在 进行 风电 功率 预 
报时 ,空气 密度 默认 为 当地 观测 的 多 年 平均 值 ,不 
再 单独 考虑 空气 密度 变化 对 风电 功率 预报 结果 的 
影响 。 风 向 也 是 影响 风电 功率 预报 的 重要 因子 之 
一 , 同 平均 风速 不同 风向 下 的 风机 产 出 功率 截然 
不 同 。 该 风电 场 80 m 高度 上 盛行 南西 南 风 和 西北 


风 ( 图 9a),33 台 风机 沿 西北 一 东南 方向 分 2 列 排 布 ， 
ECMWF 的 100 m 高 度 风 速 的 主 风 癌 为 南 风 和 西北 
风 ( 图 9b) ,与 该 风电 场 的 主 风向 比较 接近 ,但 是 在 
偏 南 风向 上 有 一 些 差异 ,这 可 能 与 模式 地 形 与 实际 
地 形 存在 一 定 差异 有 关 ”。 由 于 风机 的 偏 航 控制 
系统 使 得 风机 对 风向 有 一 定 的 自 适应 性 。 因 此 ,本 
文 的 风电 功率 预报 方程 虽然 没有 考虑 风向 因子 ,但 
是 对 风机 功率 的 预报 结果 不 会 带 来 太 大 影响 。 整 
体 上 讲 ,本 文 的 风电 功率 预报 方法 是 可 行 的 ,下 一 
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E PE: 基于 最 优 训练 期 的 PP 与 MOS 的 风电 功率 趋势 预报 对 比 


(a) 测 风 塔 80 m 高 度 风向 频率 


SSW 和 SSE 
(b) ECMWF 模 式 100 m 高 度 风向 频率 


图 9 测 风 塔 80 m 和 ECMWF 模 式 100 m 高 度 上 的 风向 频率 玫瑰 图 (单位 :%) 
Fig.9 Rose diagrams of wind directions frequency from observation at 80 m of the wind turbines 
and forecast at 100 m of the ECMWF mode (Unit: %) 


步 可 以 通过 采用 更 高 时 空 分 辨 率 的 风速 预报 产品 ， 
制作 更 精细 的 风电 功率 预报 。 同 时 ,该 风电 功率 预 
报 方 法 还 可 以 通过 引入 风向 和 密度 等 因子 来 不 断 


完善 。 
6 结论 


本 文 利 用 陕西 省 某 一 风电 场 区 域内 的 风机 观 
测 资料 以 及 ECMWF 模式 的 100m 高 度 风 速 预报 产 
品 ,基于 一 元 线性 回归 方法 ,引入 最 优 训练 期 方案 ， 
先 将 模式 的 100 m 高 度 风速 预报 订正 到 风机 轮 慌 高 
度 上 ,再 利用 PP 法 和 MOS 法 进行 风电 功率 预报 对 
比 研 究 , 主 要 结论 如 下 : 

(1) 不 管 是 风速 预报 订正 还 是 风电 功率 预报 ， 
并 不 是 最 优 训练 期 日 数 越 稳 定 预 报效 果 越 好 ,以 预 
报 日 前 1 天 之 前 的 10 d 为 评估 周期 确定 的 最 优 训练 
期 日 数 虽 然 较 不 稳定 ,但 是 其 预报 误差 较 小 ,并且 
需要 的 训练 样本 少 , 有 利于 业务 化 应 用 。 

(2) ECMWF 模式 的 100 m 高 度 风 速 预报 值 与 
风机 轮 载 高 度 处 的 风速 观测 值 之 间 的 偏差 并 不 大 , 订 
正 后 的 预报 误差 进一步 减 小 ,3~228 h 预 报 的 均 方 根 
误差 和 平均 绝对 误差 分 别 从 订正 前 的 2.2~3.9 mes"! 
和 1.7~3.2 ms 降低 为 订正 后 的 1.6~2.6 ms 和 
1.2~2.1 ms ,风速 改善 百分比 为 24%~41%。 

(3) 先 对 模式 的 100 m 高 度 风速 预报 进行 订 
正 , 再 代入 PP 法 和 MOS 法 建立 的 风电 功率 预报 方 
程 ,可 以 显著 减 小 PP 法 的 预报 误差 ,但 是 并 没有 改 
进 MOS 法 的 预报 结果 。 与 利用 风速 预报 订正 产品 
进行 风电 功率 预报 的 PP 法 相 比 ,MOS 法 可 以 省 去 风 
速 预报 的 订正 环节 而 略 胜 一 筹 。 


以 预报 日 前 1 天 之 前 的 10 d 作 为 最 优 训练 期 的 
评估 周期 ,对 利用 较 少 的 历史 资料 进行 风速 或 者 风 
电功率 预报 是 有 利 的 ,但 是 该 评估 周期 仅 适 用 于 该 
风电 场 的 33 台风 机 站 。 最 优 训 练 期 及 其 评估 周期 
的 选择 ,不 仅 与 预报 对 象 有 关 ,还 受 当 地 气候 背景 
的 影响 ,应 该 基于 不 同 评估 周期 方案 的 实验 结果 ， 
客观 确定 最 优 训 练 期 及 其 评估 周期 。 另 外 ,由 于 
ECMWF 的 分 辩 率 不 太 高 ,并 且 其 风 回 预报 与 观测 
有 一 定 差 异 , 本 文 的 研究 方法 不 适用 于 复杂 地 形 ， 
如 果 将 通过 动力 预报 方法 获取 的 高 分 辩 率 模式 结 
果 作 为 研究 资料 , 则 能 很 好 地 弥补 这 一 缺陷 。 

致谢 :感谢 陕西 省 气象 科学 研究 所 类 盼 星 博士 
和 陕西 省 气象 台 张 购 庭 博士 在 本 文 修改 过 程 中 给 
予 的 帮助 。 
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Abstract: On the basis of data sets of wind speed and wind power observed from 33 wind turbines in a wind 
farm in Shaanxi Province, China and with wind speeds at a height of 100 m forecast by the ECMWF model, the 
wind speed prediction was corrected using simple linear regression and introducing the optimal training period 
and forecast models for wind power business were comparatively studied using the perfect prognostic (PP) and 
model output statistics (MOS) methods. The optimal training period was defined as a sliding cycle, which should 
be given before the forecast date. For training days varying from 10 to 100 d, taking some past forecast days as an 
evaluation period to calculate the root-mean-square error (RMSE) of the forecast, the optimal training period was 
decided according to the variation of the RMSE of the wind speed or wind power forecast with the training days. 
A prediction equation of wind power was established using the PP method according to historical observations of 
wind speed and wind power in the training sample. Then, the numerical prediction product of wind speed was 
used to replace the observation data in the prediction equation, which was applied to a wind power prediction 
business. In contrast, the MOS method directly used the wind speed of the numerical weather prediction model to 
establish the prediction equation, which was quite different from the PP method. The result showed that wind 
speed predictions after correction and wind power were improved during the optimal training period, reducing the 
prediction error to the greatest extent. The optimal training days were determined by the RMSE of the forecast 
during the evaluation period of 10 d before the forecast date. The wind speed at a height of 100 m forecast by the 
ECMWE model was closer to that of the observation at the height of the fan hub, and the RMSE and average ab- 
solute error (ABE) respectively ranged from 2.5 to 4.1 m-s ` and 1.9 to 3.3 m»s '. The prediction error of the 
wind speed after correction became smaller, with the percentage of improvement in the RMSE of the wind speed 
forecast ranging from 33% to 46%, and the RMSE and ABE of the wind speed prediction respectively ranged from 
1.5 to 2.4 m-s ‘and 1.1 to 1.9 m-s '. Thus, wind speed predictions could be applied to wind power prediction 
businesses using both the PP and MOS methods. Although the wind power prediction error of the MOS method 
was nearly the same as that of the PP method when the wind speed prediction was corrected, the MOS method, 
which did not need to improve the quality of wind power prediction by correcting the wind speed prediction, was 
slightly better than the PP method in terms of complexity. The RMSE and ABE of the wind power forecasts of 
both the PP and MOS methods were 17%-25% and 11%-18%, respectively, within the forecast time of 3-72 h. 
Thus, both the PP and MOS methods could be applied to wind power prediction businesses. 

Key words: wind farm; wind power prediction; optimal training period; model output statistic (MOS) predic- 
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